SKIN EFFECT PROBLEM WITH THE DISPLACEMENT CURRENT 
IN MAXWELL PLASMA BY THE SOURCE METHOD 



Ph. 



Y.F. Alabina, A.V. Latyshev, A. A. Yushkanov 

Moscow State Regional University 
105005, Moscow, Radio st., 10 a 
e-mail: yf.alabina@gmail.co'm, avlatyshev@'mail.ru, yushkanov@inbox.ru 

O ; We found the analytical solution to the problem of the skin effect for 

^ ■ Maxwell plasma with the use of the kinetic equation, where the frequency 

^ ■ of electron collisions is constant. We use the specular reflection of electrons 

from the surface as a boundary condition. The behavior of an impedance 
near to a plasma resonance is considered. We consider limiting cases of skin 
effect. 

PACS numbers:52.35 - g; 52.20 - g; 52.25 - b 

Introduction. Skin effect is the plasma response to variable electromagnetic 

CnI ■ field, tangential to the surface P [2]. 

; First, the analytical solution of the skin effect problem at any value of 

(N . 

oo ■ anomaly parameter was found in [3] and [4j for plasma in metal. For the 

(si ■ gas plasma, the corresponding solution is considered in [5]. There has been 

§ : sub.tant,al interest to thi. problem B H B Um In 0. the behavior of 

plasma near the resonance is not considered. Also, research of impedance 
^ ■ with displacement current near to resonance isn't carried out in the previous 

works. For example, the displacement current isn't taken into account in [2]. 
Research of surface impedance near to plasma resonance isn't carried out in 
the rest of works. We research the behaviour of impedance near to plasma 
resonance with the displacement current in this paper. 

In this paper we continue the development of an analytical method of 
solving boundary problem for systems of equations for the electric field in 
half-space gas plasma. The basis of the method is the idea of the symmetrical 
continuation of the electric field to the conjugate half-space. We provide the 
analytical solution of the boundary problem of the skin effect theory for 



electron plasma, that fills the half-space. We formulate analytical expressions 
for electric field, distribution function of electrons and impedance. 

We assume that electromagnetic wave is incident normally to the interface 
of the plasma. In such configuration, the electric field of electromagnetic wave 
has only tangential component. We use the specular electron reflection from 
interface as boundary condition. The interface of ions on the conductivity of 
plasma is not considered. 

1. Problem statement and basic equations. 

Let Maxwell plasma fills the half-space x > 0, where x is the coordinate 
orthogonal to plasma boundary. Let the external electric field has only y 
component. Then the self-concordance electric field inside in plasma also has 
only y component Ey{x^t) = E{x)e^'^^^. Let us take the kinetic equation for 
distribution function of electrons: 

| + v.| + ei5We-g = .(/o-/(t...v)). (1) 

In (1) p is the frequency of electron collisions with ions, Cq is the charge 
of electron, /o(v) is the equilibrium Maxwell distribution function, 

/o(v) = nf^V^'exp(-/?V), ^ 



Here m is the mass of electron, Ub is the Boltzmann constant, T is the 
temperature of plasma, v is the velocity of the electron, n is the concentration 
of electrons, c is the speed of light. 

The electric field E{x) satisfies Poisson's equation 

E'\x) + ^E(x) = j yj{t, v) dh. (2) 

Assume that the intensity of an electric field is such that linear approximation 
is valid. Then distribution function can be represented in the form: 

/ = /o (1 + Cyexp(-ia;t)/i(x,//)) , 

where C = is the dimensionless velocity of electron, /j, = Cx- Let I = vtt 
be the mean free path of electrons, vt = l/v^? t = l/v. We introduce the 



dimensionless values: 

= ut, xi = -, e{xi) = — ^==E{xi). 

Below, instead of Xi we shall write again x. In new variables, the kinetic 
equation (1) and the equation on a field with the displacement current (2) 
become 

dh 

Ai— + zq h{x, ji) = e(x), zo = 1 - icur, (3) 

00 

e"{x) + Q^e{x) = -i^ [ exp(-/) hix, //') dfi', Q = -, (4) 

— oo 

(P'Ti 

where 5 = is the classical depth of the skin layer, ctq = , 

zTTuao mv 

2/2 

a = -pr, what a as the anomaly parameter. 

We formulate the boundary conditions for the distribution function of the 
electron in case of the specular electron refiection from the surface: 

/i(0,//) = /i(0, -//), 0<//<+oo. (5) 
We use the condition that function vanishes far from the surface: 

/i(+oo, //) = 0, — oo < < +00, (6) 
and conditions for electric field on the interface and far from it: 

e'(0) = e/, e(+oo) = 0, (7) 

where e/ is the given value of gradient of electric field on the plasma interface. 

So, the skin effect problem is formulated completely. We seek solution of 
system of the equations (3) and (4) in this problem that satisfy boundary 
conditions (5)-(7). 

2. The analytical solution of the problem. As a first step in the 
source method, we extend the electric field and distribution function to the 
"negative" half-space x < 0: 

e{x) = e{—x), h{x,iJ,) = h{—x, — //). (8) 

After we substitute x = to equation (8), we obtain that electric field 
and distribution function of electrons are continuous and the derivative of 



an electric field has discontinuity: e'(+0) — e'(— 0) = 2e/. In account of 
this circumstance, we introduce term with Dirac delta function to the field 
equation |10) : 



oo 

e"{x) + Q'^e{x) — 2e'g5{x) = —i^= I exp(— /i'^) h{x^ d^', 

VTT J 



(9) 



where 5{x) is the Dirac delta function. 

The third term in the left hand side of the equation (9) corresponds to 
discontinuity of a derivative of the electric field for a; = 0. 

The solution of the (3), (9), (5)-(7) can be sought as Fourier integrals (by 
variable x): 



oo 

1 



e{x) = / e'^''E{k)dk, (10) 



— oo 
oo 



h{x,ii) = ^ [ e'^''^{k,ii)dk, (11) 
27r ./ 



-oo 

oo 



5{x) = — I e'^'^dk. (12) 
27r J 

— oo 

We substitute (10)-(12) into (3) and (9). We get the following system of the 
characteristic equations: 



00 



{q'^ -k^^E{k)-2e'{<d) = -i^ J exp(-^2)$(A;,/i) dfi, 

— oo 

$(A;, fi){zo + ikjj.) = E{k). 

From these equations, we get spectral densities of the distribution function 
and electric field, respectively: 



where 



oo 

. a f exp(— /i^) dfi 



k'^ /c^i/tt J ikji + zq 



-oo 



To find the profile of the electric field in the half-space, we substitute (14) 
into (10): 



oo 



To get the function of the electron distribution in the half-space we substitute 
14) into (13). It is obvious that this spectral density is: 

Now we substitute (16) into (11). We get 

./ TOO ^^kx^^ 



h{x, fi) = — - / 



oo 



[zq + ikjj.)k'^X{k) 

3. The impedance evaluation. We introduce the dimensionless decrease 
of the electric field into the depth of plasma: 

A(a,Q) = -£M, n = ur. (17) 



From (15) this dimensionless decrement is 

ex 

A(a,n) = i J 



oo oo 

dk 2 f dk 



k^X{k) ttJ k^X{k)' 

-oo 

The value of the impedance can be calculated with the help of the formula 
see P): 

4:7riu e(x) 

Zj = 



(? de{x] 



dx 

It should be noted that 



de{x') de{x) dx ^de{x) 

Thus, according to the previous formula the impedance is 

Amojl e(0) 



Z 



c2 



or, taking into consideration the equality (17). 



Z = -i^^K{a,n). (18) 



C2 



Let us introduce parameter R, which is equal to the modulus of impedance 
in normal skin effect (when a <C 1, <C 1) 



R = 



'47TUJ 

cVo' 



where ctq is the static electrical conductivity of plasma (for a; = 0). Now 
formula (18) for the impedance can be written in the form Z = RZq, where 
Zq is the dimensionless part of impedance, 

Zo = -iy^A{a,n). (19) 

4. The analysis of the solution. We substitute the variable k = 1/t 
to the integral (19). In this case the dimensionless decrement is: 



00 

dt 



Here 



A(i/t) = 1 - - ^ y 



x{i/ty 





oo 



fi — izot 



-oo 



We will present special partial cases of formula (19). 

We will start from the case of the normal skin effect. In this case 

We suppose that there is no current of the bias. In this case the dimensionless 
impedance depends on a is 



00 

Zq = -iyJoL ^ J 



dt 



„ 1 ^ exp(— w ) du 
1 — at'^—= / ^ — 

7T U — it 



or, if we introduce 



00 



1 f exp(— w ) du 



u — iz 

oo 



the expression for the dimensionless impedance will be follows: 

oo 

dt 



Zq = -iy/a - I - 

IT J 1 



at^ to{it) 





It should be noted that for the large t: to{it) « i/t. Therefore 

00 

f dt ^ ^ ^ 1 l-i 

^0 = —lyOL — / : — — = zvQf • Res 



TT / 1 — ioii'^ z=-^ 1 — ia a/2 

— oo 

1 — i 

The expression Zq = ^ is a well-known classical result 

We will consider the anomalous skin effect in the low-frequency limit, that 
is, when a^l, r2<Cl, zq ^ 1. In this case for small t we have: 

CO 

,. , 1 f exp(—u'^)du . I— / /■ x9x ■ 1— 

toUt) « —= / « ly/TT exp( — Ut) ) « ly/TT. 

'tt u — it 



Thus the non dimensional impedance is 

00 

7 ■ ^2 /■ dt 

Zq = —i\/a- 



TT J 1 — iay^t^ 


We substitute the variable in this integral t = 1/k and obtain 

00 

.2 r kdk 



Zq = —%\fa- 



7T J — ia^Jli 



We substitute one more variable k = yoy^x, so 



i—^x/a^pK f X dx 
Zq = —i\Ja- 



TT x^ — i 



We consider the integral 

00 00 00 

J = : = Ji + iJ2, di = ~~r~r^ ^2 = 



x^ — i ' J + 1 ' J + 1 



We will calculate the integral Ji with the use of residues. The function 

(tTT \ 
— (1 + 2k) j , k = 

0,1,2. So 



oo 2 



1 f X'^dx _ 2^ TTZ/l 1 l\7r 



di = t: — 7 = TTz > Res — ; — - = — 
2 J x^ + 1 ^ z=zk 2^ + 1 6 V 



-oo 



The second integral is calculated by the decomposition of the function under 

TT 

integral into elementary fractions. As a result we have: J2 = 3^^' "^^^^5 

J = + and the expression tor the non dimensional impedance is 

3v 3 

^» = IS*' - 

which also coincides with the classical result [S]. 

Further, to study the impedance near to a plasma resonance, it is convenient 
to use the following dimensionless parameters: 

7 = — , £ = — , where (u;„ = . 

Up Up m 

Here cUp is the plasma frequency. 

We express the parameters of the problem a, il, Q through 7, e and Vc = 
vt/c, where vt = is the heat velocity of the electrons. We obtain that 

a = ^, Q = , n=-. 

e-^ 6 6 

We analyze numerically the growth of value of the modulus of impedance, 
the real, imaginary parts of impedance and argument of impedance depending 
on change of value 7 from 0.5 to 1.2 with various values of other parameters. 
If we change of the anomaly parameter a in the specified limits, the value 
7 becomes 7=1, that \s lo = Up i.e. the oscillation frequency of external 
field is the value of plasma frequency. This is plasma resonance. It would 
be interesting to consider parameters of self-concordant field near plasma 
resonance. In conclusion we show the results of numerical analysis. 

Conclusion. The analysis of plots in figure la shows that at the same 
temperature of plasma, the maximum of the modulus of impedance is reached 
at 7 = 1, i.e. for Thus, the less is the effective frequency of electrons 

collisions with particles of plasma, the greater is the modulus of impedance. 

From figure lb we see that at the same frequency of collisions of electrons 
the maximum of the modulus of impedance is reached at a; = cjp, independent 
of the temperature. Actual curves in this figure are computed for various 
values of parameter Vc = "^r/c, which depends on temperature (it is proportional 



to root square of temperature): Vc = \j2kBT j (9"m. Thus the size of the 
modulus of impedance decreases quickly with growth of temperature. 

It is interesting to note that for decrease reduction of value e from 10~^ to 
10""* (by two orders), the modulus of impedance also increases by two orders, 
more precisely 95 times. The temperature of the plasma here is 3000-fC. If 
temperature of the plasma is SOOOiT, the modulus of impedance increases by 
97 times for the same reduction in the value of s. 

For an increase in the temperature of the plasma from 3000/^ [vc = 10~^) 
to 5000-fr {vc = 13 • 10~^) the value of the modulus of impedance increases 
by 170 times, and for a change in the temperature of the plasma from 1000-fC 
to 3000-ff, the value of the modulus of impedance changes only by a factor of 
2.8. Thus, the growth of the modulus of impedance depends on the nonlinear 
change in temperature. 

In figures 2a and 2b we show the plots of the real part of impedance 
(more precisely, the plots of values Re {—Zq}). The analysis of plots in figures 
2 shows that as the frequency of collisions of electrons increases, the value 
Re(— ^o) grows at a constant temperature. If the frequency of collisions of 
electrons is constant, then this value grows as the temperature increase. 

Let us note, that in figures 1 and 2 we used the logarithmic scale for 
vertical axes. 

The analysis of dependence of argument of impedance on parameter 7 
shows that near to plasma resonance, the argument of impedance has step 
irrespective of the frequency of electron collisions (fig. 3) and from temperature 

So, the numerical analysis of plots (fig. 1, 2) shows that near to plasma 
resonance, the modulus and imaginary part of impedance have the sharp 
maximum, which is absent in low-frequency limit, or in the normal skin effect 
theory, where the argument near to resonance has step, and the real part of 
impedance has the sharp maximum. 




Fig. la. The modulus of impedance. For 
tlie curves 1,2,3, e = 10"^ 10"^, lO^^ 
respectively; and Vc = 10~^. 



Fig. lb. The modulus of impedance. For 
the curves 1,2,3, v,, = 6 ■ 10~^ lO'^, 13 • 
10~^ respectively; and e = 10~^. 
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Fig. 2a. The real part of impedance. For 
the curves 1,2,3, e = 10"^ 10"^ 10"^ 
respectively; and Vc = 10~^. 
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Fig. 2b. The real part of impedance. For 
the curves 1,2,3, = 13 ■ IQ-^, 10"^ 6 • 
10~^ respectively; and e = 10~^. 
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Fig. 3. The argument of impedance. For the curves 1,2,3, £ = 10 ^, 10 ^, 10 ^ 

respectively; and Vc — 10~^. 
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